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ABSTRACT 

We present the results of a comparative study of HI narrow self-absorption 
(HINSA), OH, ^^CO, and C^^O in five dark clouds. We find that the HINSA 
generally follows the distribution of the emission of the carbon monoxide iso- 
topologues, and has a characteristic size close to that of ^^CO. This confirms 
earlier work (Li & Goldsmith 2003) which determined that the HINSA is pro- 
duced by cold HI which is well mixed with molecular gas in well-shielded regions. 
The OH and ^^CO column densities are essentially uncorrelated for the sources 
other than L1544. Our observations indicate that the central number densities of 
HI are between 2 and 6 cm~^, and that the ratio of the hydrogen density to total 
proton density for these sources is 5 to 27x10"^. Using cloud temperatures and 
the density of atomic hydrogen, we set an upper limit to the cosmic ray ioniza- 
tion rate of 10^^^ s^^. We present an idealized model for HI to H2 conversion in 
well-shielded regions, which includes cosmic ray destruction of H2 and formation 
of this species on grain surfaces. We include the effect of a distribution of grain 
sizes, and find that for a MRN distribution, the rate of H2 formation is increased 
by a factor of 3.4. 

Comparison of observed and modeled fractional HI abundances indicates ages 
for these clouds, defined as the time since initiation of HI H2 conversion, to 
be 10^'^ to 10'^ yr. Several effects may make this time a lower limit, but the 
low values of ylhi we have determined make it certain that the time scale for 
evolution from a possibly less dense atomic phase to almost entirely molecular 
phase, must be a minimum of several million years. This clearly sets a lower 
limit to the overall time scale for the process of star formation and the lifetime 
of molecular clouds. 

Subject headings: ISM: atoms - individual (hydrogen) 
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1. INTRODUCTION 

The study of atomic hydrogen in dense interstellar clouds has a long history. Very 
early measurements of dark clouds indicated no correlation of atomic hydrogen and dust 
opacity (Bok, Lawrence, & Menon 1955). Subsequent measurements in fact showed an 
anticorrelation between the HI column density and extinction (e.g. Garzoli & Varsavsky 
1966). These pioneering observations of the Taurus region were followed by a broader survey 
of dust clouds by Heiles (1969) in which a number of spectra showing "fluctuations" were 
found, but one object showed clear evidence for "self-absorption": a narrow absorption 
feature superimposed on a wider emission line. In this source in Taurus, the velocity of the 
HI self-absorption agreed with that of the OH emission. 

Since that time, a number of large-scale and focused studies of HI in dark clouds have 
been carried out (e.g. Knapp 1974; Heiles & Gordon 1975; Wilson & Minn 1977; Mc- 
Cutcheon, Shuter, & Booth 1978; Myers et al. 1978; Bowers, Kerr, & Hawarden 1980; Minn 
1981; Batrla, Wilson, & Rahe 1981; Shuter, Dickman, & Klatt 1987; van der Werf, Goss, 
& Vanden Bout 1988; Feldt & Wendker 1993; Montgomery, Bates, & Davies 1995). With 
improved sensitivity and higher angular and spectral resolution, most of these studies iden- 
tify self-absorption features. The inherent ambiguity between a temperature fluctuation in 
the emitting region and a distinct, colder, foreground source responsible for the absorption 
makes it difficult to distinguish between these two quite different situations for an arbitrary 
line of sight (e.g. Gibson, Taylor, Higgs, & Dewdney 2000; Kavars et al. 2003), and it is 
still only for relatively nearby clouds that the location of the cold atomic hydrogen can be 
identified with certainty. 

A key element in pinpointing the location of cold HI is its association with molecular 
emission, and the greatest difference between early and current observations of HI self- 
absorption is the recognition that optically-identified dust clouds are cold (T ~ 10 K), dense 
(n > 10"^ cm~^) primarily molecular regions. The availability of reasonably well-understood 
tracers which are surrogates for the largely-unobservable H2 means that determining the HI 
to H2 abundance ratio in nearby dark clouds is now practical. These regions are of interest in 
their own right, as they contain the dense cores in which low mass star formation can occur, 
and in regions such as the Taurus molecular cloud, these are the only regions in which star 
formation appears to be occurring. Additionally, because the visual extinction of these dark 
clouds is sufficient to make photodestruction of H2 unimportant in their central regions, they 
are ideal laboratories for comparative studies of the atomic and molecular forms of hydrogen 
under relatively simple conditions. 

A recent survey of dark clouds using the Arecibo telescope revealed that over 80% of 
the clouds in Taurus exhibit HI narrow self-absorption (Li & Goldsmith 2003), which we 
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gave the acronym HINSA. The atomic hydrogen responsible was shown to be associated 
with molecular material traced by OH, C^^O and ^^CO on the basis of (1) the agreement in 
velocities, (2) the close agreement of nonthermal line widths, and (3) the low temperature 
implied by the depth of the HI features and the small values of the nonthermal contribution 
to the line width. Mapping data for only one cloud was available, which showed reasonable 
agreement between the column densities of cold HI, and OH, and the integrated intensity of 
Ci«0. 

Wc have continued and extended this work, and in the present paper report observations 
of four clouds: L1544, B227, L1574, and CB45. In Section 2 we present the observations of 
HI absorption, OH, C^'^O ^'^CO and CI emission. In Section 3, we discuss the spectra and 
determination of the column densities of the species we have observed. We also present the 
maps of these tracers in each of the sources. In Section 4 we present and discuss correlations 
of the column densities of different tracers. In Section 5 we derive the number densities of 
the species and determine the HI to H2 ratio for each of the clouds. In Section 6 we present 
models for the time dependence of the atomic and molecular hydrogen density. Comparison 
of models and our data allows us to draw conclusions about the lifetime of these clouds in 
Section 7. 



2. OBSERVATIONS 

The HI and OH data were obtained using the 305 m radio telescope of the Arecibo 
Observatory^. The data were obtained in four observing sessions between November 2002 
and May 2003. The telescope was calibrated as described by Li & Goldsmith (2003) and 
both the HI and OH data, which were obtained simultaneously using the "L-widc" receiver, 
have been corrected for the main beam efficiency of 0.60. The elliptical beam has FWHM 
beam dimensions of 3'.1 by 3'. 6. We observed using total power "on source" observations 
only, and removed the receiver noise and system bandpass by fitting a low-order polynomial 
to the spectrometer output. For these observations, the autocorrelation spectrometer was 
configured to have three subcorrelators with 1024 lags each covering 0.391 MHz centered on 
frequencies of the HI, OH(1665 MHz) and OH(1667 MHz) transitions. A fourth subcorrelator 
with 1024 lags covering a 0.781 MHz bandwidth was centered on the HI line, and was used 
for removing the instrumental baseline when wider emission features were encountered. The 
nominal channel spacing for the HI observations was 0.08 km s~^, but data were typically 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated 
by Cornell University under a cooperative agreement with the National Science Foundation. 
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smoothed to a velocity resolution of 0.16 km s ^. 

The receiver temperature was typically 30 K at 1420, 1665, and 1667 MHz, but the total 
system temperature was more than double this value when observing 21-cm emission lines. 
The integration time per position was 1 minute. This was sufficient to give good signal to 
noise for the HI absorption features, and also satisfactory OH 1667 MHz emission spectra. 

Observations of ^^CO, ^^CO, and C^^O were obtained using the Five College Radio 
Astronomy Observatory 13.7 m radio telescope ^, during April and May 2003 and January 
2004. The 32 pixel Sequoia receiver was configured to observe ^^CO and C^^O simultaneously 
for most of the data taking. The autocorrelation spectrometer for each isotopologue had a 
bandwidth of 25 MHz and 1024 lags, corresponding to a channel spacing in velocity of 0.067 
km s~^at 110 GHz. We mapped each source using "on the fly" data taking, covering a region 
determined by the HI observations obtained previously. The angular size of the maps ranged 
from 12' to 27'. System temperatures were typically 200 K for ^^CO and C^^O and 450 
K for ^^CO. The integration time per sample point was nominally 10 s. The oversampled 
data was convolved to a circular Gaussian having the same solid angle as the Arecibo beam, 
corresponding to a FWHM equal to 3.3 '. We smoothed the data in velocity to the same 0.16 
km s~^ resolution used for the HI and OH. Additional data on ^^CO for the source L1544 
were obtained to be able to make accurate corrections for the optical depth of the ^^GO 
emission. All FCRAO data were corrected for appropriate main beam efficiencies, which are 
0.49 for i^CO and G^^O and 0.45 for ^^CO. 

We obtained the data on atomic carbon (GI) using the Submillimeter Wave Astronomy 
Satellite (SWAS; Mclnick et al. 2000). Observations having total integration time of 1 hour 
per position were obtained using position switching to a reference position selected to have 
very weak or undetectable ^^GO J^O emission. The data have been corrected for the main 
beam efficiency of 0.90. The SWAS beam size is relatively large (3'. 5 x 5'.0 at the 492 GHz 
frequency of the GI line) and combined with low intensities of this transition in dark clouds, 
we were only partially successful in being able to delineate the distribution of GI relative to 
that of the other species studied. These results are further discussed in Section 3. 



^The Five College Radio Astronomy Observatory is operated by the University of Massachusetts with 
support from the National Science Foundation. 
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3. COLUMN DENSITIES, SPECTRA, AND DISTRIBUTIONS 

Determination of the distances to dark clouds such as those in the present study presents 
considerable difficulties. Lf544 is associated with the Taurus dark cloud, for which the dis- 
tance given here is the average of various techniques described by Elias (1978). B227 is at 
a galactic latitude of -0.46", and is located in a complex region of considerable extinction. 
Arquilla (1984) observed the same cloud as studied here, and with a B-V color excess tech- 
nique, derived a distance of 600 pc but there is obviously a large uncertainty. Tomita, Saito 
& Ohtani (1979) studied the dark nebula L1570 which is often (e.g. Clemens & Barvainis 
1988) used interchangeably with B227, but there is a difference of over 8' between the Ar- 
quilla and Tomita et al. positions. The latter authors, using star counts derived a distance 
300 pc, but again with significant uncertainty. Earlier, Bok & McCarthy (1974) gave an 
"adopted" distance of 400 pc for B227. We here adopt a distance of 400 pc. L1574 and CB45 
are lumped together by Kawamura et al. (1998) and are assigned the 300 pc distance of the 
not too distant L1570 determined by Tomita, Saito, & Ohtani (1979), which we adopt here. 
The source coordinates and distances we have adopted are given in Table 1. 



3.1. Cold HI 

The column density of the cold, relatively quiescent atomic hydrogen was found by as- 
suming an excitation temperature of 10 K, and solving for the optical depth of the absorbing 
gas as described in Li & Goldsmith (2003). The maximum values of Tcouhi are modest, 
ranging from 0.34 for L1574r to 0.69 for L1574b. As discussed in the earher paper, the effect 
of foreground HI is primarily to reduce the apparent value of the 21-cm absorption optical 
depth. We have corrected the optical depth as described in Li & Goldsmith (2003), with 
values of p determined from the distances given in Table 1. The correction was less than 
15 %. The corrected optical depth and the line width of the gaussian fit to the absorption 
profile together yield the column density of cold absorbing gas using equation 18 of Li & 
Goldsmith (2003). 

We present the peak main beam temperature Tmf, at the positions of maximum emission 
of each cloud in Table 2. The values given are the peak antenna temperatures divided by 
the main beam efficiency. 



-6- 



3.2. OH 

The OH emission is almost certainly optically thin. The peak main beam temperatures 
are typically 0.5 K, except in L1544 where the line is a factor of three stronger. The modest 
A-coefficient (7.8xl0~^^ s~^; Destombes et al. 1977) and typical coUisional rate coefficients 
(> 2xl0~^^ cm^ s~^; Dewangan, Flower, & Alexander 1987; Offer & van Dishoeck 1994) 
result in a critical density < 4 cm~^, so that the excitation temperature even in these cold 
clouds will be close to the kinetic temperature. We use equation 15 of Li & Goldsmith 
(2003) to calculate the OH column densities from the integrated intensities. 

3.3. Carbon Monoxide Isotopologues 

The C^^O 1-0 line is relatively easy to excite and the excitation temperature of this line 
is again expected to be close to the kinetic temperature. In our sample, the C^^O antenna 
temperatures are sufficiently small compared to the ~ 10 K kinetic temperatures to justify 
the assumption of low optical depth. We employ equation 18 of Li & Goldsmith (2003) to 
derive the total column density, assuming LTE and an excitation temperature of 10 K for 
all sources except L1544 (disucssed below). The ^^CO emission is also optically thin in all 
sources except LI 544, where the positions of stronger emission near the center of the cloud 
have optical depth of order unity. The comparison of ^^CO and C^^O emission (discussed 
further below) is consistent with the generally optically thick ^^CO except in the central 
region of LI 544. 

In order to assess and correct the situation for this cloud, we obtained a map of the 
^^CO emission from L1544. We derived the excitation temperature from the main beam 
temperature, assuming that the emission is optically thick. Values for this source ranged from 
9.5 K to 13 K, a few K warmer than typical for dark clouds. These excitation temperatures 
were used to compute the C^^O column densities. We assumed that the ^^CO has the same 
excitation temperature, and from its main beam temperature, we derived the actual optical 
depth and corresponding column density. The maximum ratio between the actual and the 
"thin" optical depths is a factor of 1.75, so the effect of the saturation on the ^^CO is 
relatively modest and confirms the appropriateness of assuming low optical depth for the 
C^^O and for the ^^CO in clouds other than L1544. 

In Table 3 we give the column densities for the cold HI, OH, ^^CO, and C^^O at the peak 
HINSA positions of the clouds observed. We also give the column density of H2 obtained 
from that of the carbon monoxide isotopologues using the expressions of Frerking, Langer, & 
Wilson (1982) for Taurus. The molecular hydrogen column densities for each cloud derived 
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using ^^CO and C^^O agree quite well. 

3.4. Uncertainties in Determination of the Column Densities 

The dominant contributions to the uncertainties in the integrated intensities of the 
various spectral hues observed in this study are quite different. For the cold HI, this issue 
is discussed at length in Li & Goldsmith (2003); for the HINSA, there are two major 
contributors. The first is the foreground correction, which requires knowledge of the local HI 
distribution and the distance to the cloud for evaluation. This correction contributes at most 
15% for the present sample of clouds. The second contributor is the fitting of the background 
spectrum, which adds about 10% for narrow absorption lines. The derived column density 
is inversely proportional to the assumed gas kinetic temperature. The arguments in Li & 
Goldsmith (2003) indicate that in these clouds, the HI responsible for the absorption is in 
cold gas traced by the carbon monoxide. The kinetic temperature is assumed to be 10 K, 
but the uncertainty is about 30%. The total uncertainty for HINSA is thus on the order of 
50%, with the statistical errors making a negligible contribution. 

For the much weaker OH lines, statistical errors are more important, especially in the 
outer portions of the clouds. The column density is only very weakly dependent on the 
kinetic temperature, and the emission is almost certainly optically thin. As discussed in 
(Li & Goldsmith 2003) ignoring the background can produce a systematic underestimate 
of the OH column density by a factor ~ 2, which has not been corrected here. For the 
carbon monoxide isotopologues, the uncertainty in the kinetic temperature dominates the 
determination of the column density (again except for the cloud edges). Assuming LTE, 
the column density changes by approximately 15% if the kinetic temperature varies from a 
nominal value of 10 K to 7 K or 15 K. As discussed above, for the central positions of L1544 
the ^^COis shghtly optically thick. The correction, which has been included, depends itself 
on the kinetic temperature, and we estimate that this introduces an additional uncertainty 
of 25%, but only for this source. We estimate that there is a factor of approximately 3 
uncertainty in the H2 column densities due to variations in the fractional abundance of the 
carbon monoxide isotopologues and uncorrected excitation and radiative transfer effects. 

3.5. Spectra 

In Figure 1 we show the spectra of four species at the reference position for each of 
the four clouds we have studied. As discussed in detail by Li & Goldsmith (2003), the 
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C^^O and ^^CO lines agree well in velocity and nonthermal line width with the OH emission 
and HINSA features. We here see, particularly evident for CB45, the close agreement of 
the complex absorption profile of the HINSA with the profiles of the three emission lines 
observed in the direction of this cloud. This strongly supports the contention that the cold 
HI is sampling the same region of space as the emission from the species studied here. The 
HI spectrum in the direction of B227 is complex, with various peaks and dips. Only the one 
of these, which is the narrowest and most clearly defined absorption feature, is associated 
in velocity with detectable molecular emission. This illustrates the difficulty of determining 
through HI spectroscopy alone whether any particular HI feature is a result of a variation in 
the temperature of the emitting region or, if a dip, is due to absorption by cold gas in the 
foreground. 

In the direction of L1574 there are two narrow HI absorption features; each is quite 
well locahzed and has a very distinct morphology. The velocity of the "blue" feature shown 
here at the (0,9) ^ position, the location of its maximum emission intensity, is km s~^, 
while the velocity of the "red" feature, which peaks at the (0,0) position, is 3.5 km s^^. In 
what follows we will refer to the cloud defined by the km s^^ feature as L1574b and that 
defined by the 3.5 km s~^ feature as L1574r. While the emission from L1574b overlaps that 
of L1574r, there is no clear indication of any physical connection, so we consider these to be 
two independent clouds at essentially the same distance along the hue of sight. 



^All offsets in this paper are given as (minutes of arc in RA, minutes of arc in Decl.) 




Fig. 1. — Spectra of different species in the four clouds included in this study. L1544 and 
B227 spectra are at the nominal central position of each cloud while the offsets for LI 574 
and CB45 are in minutes of arc offset relative to the central position. The offset position for 
L1574 corresponds to the peak of the km s~^ feature and is denoted L1574b. The spectra 
in for each position are arranged such that the HI hnes, having the greatest intensity, are 
plotted after scaling by factor 1.0. The ^^CO hnes are plotted after scahng by a factor of 5 
for L1544, 20 for B227, 25 for L1574, and 10 for CB45, and are the strongest after those of 
HI . The OH(1667) spectra are plotted as dotted lines after scaling by factors 15 for L1544, 
50 for B227, 25 for L1574, and 40 for CB45, and fah below those of ^^CO. The C^^O Hnes 
are plotted after scahng by factors 5 for L1544, 20 for B227, 25 for L1574, and 10 for CB45, 
and are the weakest spectral features. 
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3.6. Maps of Distributions 

In Figures 2 through 5 we show the HINSA cold-HI column density (color) and ^'^CO 
and C^^O column densities (contours) in the four clouds. The millimeter data have been 
smoothed to approximate the spatial resolution of the Arecibo telescope, as described pre- 
viously. Comparison of the distributions of the molecular emission, particularly the carbon 
monoxide isotopologues, reinforces the conclusion of Li & Goldsmith (2003) that the HI 
producing the narrow self-absorption features is in the cold, well-shielded interiors of these 
molecular clouds. 

In general the morphology of the molecular emission and HI absorption are quite similar. 
This is particularly evident if one compares the atomic hydrogen with the ^^CO in L1544 and 
L1574. In the latter cloud, we see very good agreement in both of the velocity components. 
The emission from CB45 peaks far from the (0,0) position, which is taken from Clemens 
& Barvainis (1988). Systematic mapping around the nominal central position revealed 
the local maximum found within a few arcminutes of the (0,0) position, but indicated an 
extension to the north. Following this up led to the far stronger maximum centered at the 
(-3,18) position, which is close to the center of the HINSA absorption, although we see in 
Figure 5 that the HI is concentrated in two features located on either side of the C^^O 
emission maximum. 

The exact peak of the C^^O is offset from the strongest HINSA feature in two of the 
four sources, but only by a small fraction of the size of the cloud. We thus feel that (1) the 
general distributions are very similar; (2) the strongest molecular emission coincides quite 
closely with the strongest absorption by cold HI; and (3) there is no indication of any limb 
brightening or other suggestion that the cold HI producing the HINSA is not well mixed with 
the molecular gas. As discussed by Li & Goldsmith (2003), this may be due in part to the 
reduced absorption coefficient of HI as the temperature increases, and thus we cannot exclude 
the possibility that there is significant warm HI in an envelope or photon dominated region 
(PDR) surrounding the molecular cloud. However, based on comparison of line profiles and 
distributions of line intensities, the narrow line HI absorption we are studying here is best 
interpreted as being produced by cold HI which is well mixed with the molecular gas. 
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L1544 




Offset (arcminutes) 

Fig. 2. — Maps of the cold HI column density producing narrow absorption features (HINSA 
in color; scale on right), contours of column densities of ^^CO (blue), and of C^^O (magenta) 
in L1544. The contours are at 30, 50, 70, and 90% of maximum values, which are 7.0 x 10^^ 
cm-2 for and 1.2 x 10^^ cm-^ for C^^O. 
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B227 
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Fig. 3. — HINSA, ^^CO and C^^O column densities for B227. The color scale and relative 
contours are the same as in Figure 2. The maximum value of the column density is 3.1 x 10^^ 
cm-2 for i^CO and 0.26 x 10^^ cm'^ for C^^O. 
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L1574 Red L1574 Blue 
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Fig. 4. — HINSA, ^^CO and C^^O column densities for L1574. The color scale and relative 
contours are the same as in Figure 2. For L1574b (velocity = km s~^), the maximum 
values of the column density are 3.6 x 10^^ cm~^ for ^^CO and 0.51 x 10^^ cm~^ for C^^O. 
For L1574r (velocity = 3.5 km s~^), the maximum ^^CO column density is 0.94 x 10^^ cm"^. 
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CB45 




-5 -10 

Offset (arcminutes) 



Fig. 5. — HINSA, ^^CO and C^^O column densities for CB45. The color scale and contours 
are the same as in Figure 2. The maximum values of the contours are 6.2 x 10^^ cm~^ for 
i^CO and 0.68 x 10^^ cm'^ for C^^O. 
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3.7. CI 

The only source sufficiently extended for its distribution to be probed by SWAS is 
Li544. In Figure 6 we show the integrated intensity of CI and the column densities of 
^^CO and cold HI (HINSA) along a cut in Right Ascension through the center of L1544. As 
indicated by the two-dimensional maps presented above, the ^^CO and the cold HI basically 
track each other. The CI, however, shows a much less centrally-peaked distribution. Its 
intensity has dropped by a factor of less than two from the center to the edge of the map, 
while the ^^CO and cold HI have diminished by a factor ~ 5. While the optical depth 
of CI may not be negligible, this data is strongly suggestive that the CI is predominantly 
in a relatively extended structure, and that the dense core with associated large column 
density seen in most molecular species does not dominate the CI emission from this dark 
cloud. This picture is supported by the relatively large line widths of the CI. As discussed 
in Li & Goldsmith (2003), these are generally much broader than those of HINSA or the 
carbon monoxide isotopologues. As shown in Figure 7 of that paper, the ^^CO and HINSA 
line widths, which are largely nonthermal, are well correlated, while the CI line widths are 
essentially uncorrelated with those of HINSA and molecular species. In the case of L1544 our 
improved data here^ give a CI FWHM line width of typically 1.5 km s~^. Comparison with 
the ^ 0.5 km s^^ FWHM line widths of C^^O and ^'^COindicates that the more extended 
region in which the CI is abundant has a larger turbulent velocity dispersion than does the 
region in which the molecular emission and HINSA absorption are produced. 



^The line width given here is narrower than that given in Li & Goldsmith (2003), presumably due to the 
higher signal to noise ratio. 
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Fig. 6. — Cut through the center of L1544 showing column density of ^^CO (divided by a 
factor Ix 10^^), the column density of cold HI producing the narrow HINSA feature (divided 
by a factor 1x10^^), and the integrated intensity of the CI ^Pi - "^Pq transition (multiplied 
by a factor of 4). The CI emission is clearly more spatially extended than the molecular 
emission or the HI absorption. 
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4. CORRELATIONS OF COLUMN DENSITIES 

4.1. Carbon Monoxide Isotopologues 

In order to trace as well as possible the molecular content of these clouds, we have com- 
pared the emission from ^^CO and C^^O. We see from Figure 7, where we plot N(^^CO) /N(C^^O) 
as a function of the C^^O column density, that this ratio decreases systematically as N(C^^O) 
increases, or as we move from the outer to the inner portions of the clouds. As mentioned 
above, the ^^CO emission in the central part of L1544 suffers from moderate saturation, and 
we have corrected the optical depth and column density for this effect, as discussed above 
in Section 3. 
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Fig. 7.— Ratio of the ^^CO to the C^^O column density as a function of N(C^^O). There is 
no data for the red component of L1574 as C^^O has not been rehably measured in the 3.5 
km s~^ component of the cloud. The points for L1544 have been corrected for saturation as 
discussed in the text. 
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It is possible that the strong dependence of N(^^CO)/N(C^^O) on the column density 
of C^^O is a result of uncorrected saturation in the more abundant species, but this seems 
unlikely for the following reasons. First, the L1544 data after correction for modest saturation 
(a factor < 2 for the most saturated positions) agrees well with that of the other sources. 
Second, the line temperatures at the central positions having the strongest emission in the 
three other sources (Table 2), are significantly below the expected kinetic and J = 1— >0 
excitation temperatures. Thus, saturation appears unlikely for the central positions of the 
three other sources, and very unlikely for the outer parts of these clouds. Third, the actual 
values of the ^^CO to C^^O ratio at the central positions of L1544, CB45, and L1574b (the 
peak of the blue velocity component at offset (0,9) relative to the cloud center) are all in 
the range 6 - 10, while the center of B227 has a ratio of 12.6. These values are basically 
consistent with the expected oxygen and carbon isotopic ratios (Langer & Penzias 1990, 
1993). 

Our interpretation is that the -"^^CO to C^^O ratio is being increased by isotope-selective 
processes operative in the outer parts of these dark clouds. Processes include chemical iso- 
topic fractionation (e.g. Watson, Anicich, & Huntress 1976) and isotope-selective photodis- 
sociation (e.g. Bally & Langer 1982). The chemical isotopic fractionation is particularly 
effective in cold regions, but requires the presence of significant fraction of carbon in ionized 
form to operate. The isotope-selective photodissociation relies on reduced rates for this 
process resulting from self-shielding in line destruction channels. Studies including both 
processes indicate that chemical isotopic fractionation will be the dominant mechanism in- 
creasing the ^^CO to ^^CO ratio and the ^^CO to C^^O ratio at visual extinctions between 
0.5 and 1 mag, if the gas temperature is ~ 15 K (Chu & Watson 1983; van Dishoeck & 
Black 1988). Both may work to increase the ^"^CO to C^^O ratio and both become more 
effective in regions of reduced extinction where either the abundance of ionized carbon or 
the rate of photodestruction of carbon monoxide isotopologues is greater. 

Plotting the ratio as a function of N(^^CO) also shows a strong inverse correlation, 
similar to those of Langer et al. (1989) for a single cloud, B5. Those authors found a 
maximum of the ^^CO to C^^O ratio of approximately 30, close to the peak in Figure 7. We 
do not, however, find the drop in the column density ratio seen at the very outer edge of the 
cloud found by Langer et al. (1989). We feel that the results presented here suggest that 
^^CO is the superior probe of column density in the outer portions of these clouds. This is 
because the fractionation enhances the abundance of ^^CO and thus compensates for the 
overall drop in the molecular abundances in the outer portions of the cloud resulting from 
photodestruction. The good correlation of ^^CO with visual extinction in the outer parts of 
clouds is further confirmation of this fortuitous cancellation of opposing effects. 
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4.2. OH 

We have obtained OH emission spectra simultaneously with those of atomic hydrogen, 
and in Figure 8 we present a comparison between the OH column density and that of 
^^CO. For the sources other than L1544 there is no clear correlation, and while there is 
some segregation of the various clouds in terms of their total ^^CO column density, there 
is no distinction to be made among them in the column density of OH. While a detailed 
explanation is beyond the scope of this paper, it does suggest that a large fraction of the 
OH emission resides in a "skin" of material in each cloud, and thus once a certain minimum 
column density (or visual extinction) is reached, the amount of OH present along the line 
of sight is essentially independent of the total column density. N(OH) shows a large scatter 
which may reflect the environment, the cloud density, or some combination of effects. 

L1544 is the single cloud which the column densities of the two species are clearly 
correlated. The column densities are related as <N(OH)> = 10^^'^ <N(^^CO)> and N(OH) 
oc [N(^'^CO)]*^'®^°''^. The difference in behavior between this cloud the others observed may 
be related to its higher density and/or the more evolved state of this source. 

4.3. Carbon Monoxide and Cold HI 

In Figure 9 we show the correlation of the column densities of ^^CO and cold HI in the 
four sources included in this study. 
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Fig. 8. — Column density of OH plotted as a function of the column density of CO. The 
^^CO data for L1544 have been corrected for saturation as discussed in the text. 
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Fig. 9.— Plot of the column density of cold HI (HINSA) versus that of ^^CO. The ^^CO 
data for L1544 have been corrected for saturation. The three diagonal lines are not fits but 
represent ratios N(HINSA)/N(^3qo) = 10^-° (highest), lO^-^, and lO^-o (lowest). 
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Although there are some outlying points, it is clear that the column density of cold HI 
is highly correlated with that of ^^CO. It is also clear that there is some difference in the 
ratio of the two column densities from one cloud to another. L1544 has the lowest ratio, 
~ 650, B227, CB45, and L1574b are fit quite well by a ratio of 3100, and L1574r, by a 
ratio of 6300. Converting the ^^CO column densities to total molecular column densities 
is complicated by the varying extinction that probably characterizes these different lines of 
sight. For reference, if we take a characteristic ratio N(^^C0)/N(H2) ~ 1.4x10"^ (the value 
for Taurus from Prerking, Langer, & Wilson 1982, omitting the offset term), the highest 
line in Figure 9 correponds to a ratio N(HI)/N(H2) = 0.014, and the lowest line to a ratio a 
factor of 10 smaller. 

Thus, while there is approximately a factor of 10 difference in N(HINSA)/N(^^CO) 
between the cloud with the smallest and that with the largest value of this ratio, there is a 
clear correlation within each cloud which covers a range of ~ 20 in the column density of 
either species, and for the ensemble of clouds, covering a range of 100 in column density. 

The values for the fractional abundance of HI that we have obtained are close to but 
somewhat larger than expected in a steady-state, and the implications of this are discussed 
in Section 7. 



5. ATOMIC AND MOLECULAR DENSITIES 

In order to extract the maximum value from the data about atomic hydrogen in molec- 
ular clouds, it is helpful to determine the actual number densities of different species. To 
do this, we model the density of HI, C^^O, and ^^CO in each cloud as a triaxial Gaussian, 
characterized by a peak value and three widths. We determine the two widths on the plane 
of the sky from the extent of the emission between half power points along two perpendicular 
axes. These are taken along (albeit imperfect) axes of symmetry for L1544, and along the 
directions of Right Ascension and Declination for the other three clouds. 

The use of a Gaussian is reasonable given the column density distributions and the 
moderate dynamic range of our data. Given the ~ 3' angular resolution, we do not see any 
strong central peaking characteristic of the power law distributions often found for cloud 
cores; the small, dense regions are presumably significantly beam diluted in the present 
study. A Gaussian function makes it relatively easy to model the significant asymmetry 
seen to be characteristic of our sources, which is difficult to do with a power law density 
distribution. Using the distances given in Table 1, we obtain the FWHM source dimensions 
given in Table 4. The FWHM along the line of sight is taken to be the geometric mean of 
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the two measured dimensions. The central density of species X is then given by 

riceniX) = — —r . (1) 

l^ZpWHM LOS\-^ ) 

The HINSA and ^^CO sizes for L1544 arc lower limits due to incomplete coverage of the 
emission. This cloud is immersed in a region of lower column density (Sncll 1981) which 
extends to the northwest. The more extended maps suggest that the half power dimensions 
in this direction are only slightly greater than our limits. We see that the cloud dimensions 
in HINSA and ^"^CO emission are comparable, which suggests that these two species are 
coextensive. This is consistent with the fact that the ^^CO nonthermal line widths are 
relatively close to those of HINSA (Li & Goldsmith 2003). 

The regions of C^^O emission are smaller than those of HINSA or ^^CO, but by a factor 
less than two. The C^^O line widths are correspondingly smaller than those of HINSA. 
Thus, while the two molecular tracers both delimit regions comparable in size to that of 
the cold HI, it appears that the ^^CO is probably the better tracer for comparison with the 
atomic hydrogen to determine fractional abundances. The central densities given in column 
5 of Table 4 are intended to be representative of the central portion of the cloud, but given 
the moderate angular resolution of our observations and the fairly gradual fall-off of the 
densities, it is characteristic of a significant fraction of the cloud volume. 

In order to determine the H2 densities in the centers of the cloud, we have to adopt 
fractional abundances for C-'^^O and ^^CO; these inevitably are somewhat uncertain. In 
particular, depletion of carbon monoxide isotopologues onto dust grains, while not so evident 
as for some other species, is a potential problem for the densest, most quiescent regions of 
dark cloud cores (e.g. CascUi et al. 1999; Bacmann et al. 2002; Bergin et al. 2002; Tafalla et al. 
2002; Redman et al. 2002). However, all of these studies indicate that significant depletion of 
rare isotopologues of carbon monoxide is restricted to the central, densest regions of clouds. 

The map of CCS emission in L1544 (Ohashi et al. 1999) shows an elliptically shaped 
region approximately 1'.5 by 3' in extent. There is evidently a central region in which is 
molecule is significantly depleted, but it is between 1' and 1'.5 in diameter. A similar-sized 
region is heavily depleted in L1498 (Willacy, Langer, & Velusamy 1998). For more distant 
clouds, regions of similar physical size (1.3 - 3 xlO^^ cm) would appear smaller. These 
"depletion zones" are heavily beam diluted (by a factor ~ 10 in the single beam looking at 
the center of the cloud), and would of course be irrelevant for the vast majority of positions 
we have studied. 

The region over which depletion has a major impact can also be seen in terms of the 
density required for significant depletion, which is 10^ - 10^ cm~^ (Bacmann et al. 2002; 
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Tafalla et al. 2002). This is basically two orders of magnitude greater than the average 
densities we find, consistent with the depletion being confined to a small portion of the line 
of sight as well as of the 3'. 3 beam. Using empirically determined density profiles, Tafalla et 
al. (2002) find that characteristic radii within which this density is attained are ~ 10^^ cm, 
or 1'. 

The depletion rate is the rate of collisions which result on sticking to a grain. For 
molecules, the sticking coefficient S is taken to be unity (Burke & Hollenbach 1983), while 
for HI (see Section 6.1 below), S = 0.3. The depiction rate is proportional to the sticking 
coefficient and inversely proportional to the square root of the particle's mass (e.g. Draine 
1985). The depletion time scale is the inverse of this rate, which gives 

^CO depletion 
'^HI depletion 

Scaling the results of Section 6.1 for molecular hydrogen formation to atomic hydrogen 
depletion and then to carbon monoxide depletion, we obtain tco depletion = 2.1xl0^/n//2 yr- 
This is consistent with the results of Caselli et al. (1999) and Bacmann et al. (2002), that 
Tco depletion — 10'^ yr for uj/j = 10^ cm~^. The relatively short time scale for molecular 
depletion is thus a result of the very high densities in the cloud cores which overall have 
shorter evolutionary time scales than those characteristic of the lower density clouds in 
which they are embedded. 

The modest effect of the depletion on determining the column density even in the 
beam including the cloud core is illustrated by the result of Caselli et al. (1999) who derive 
that the "missing mass" of H2 from depletion of tracers is 2.3 solar masses in L1544. The 
column density in the central 3'. 3 beam assuming undepleted abundance of C^^O is 4 solar 
masses. Thus, we could be underestimating the H2 column density by as much as 50% in 
this direction. This would have the effect of decreasing the HI to H2 abundance ratio in 
this (and possibly other) core directions, but would not change to an appreciable degree the 
correlation of HI and H2 or the implications of the ratios determined. 

Using the molecular hydrogen column densities and the source sizes, we derive the 
central densities of atomic and molecular hydrogen given in Table 5. We see that the atomic 
hydrogen densities arc between 2 and 6 cm^"^. The density of molecular hydrogen derived 
from C^^O is on the average a factor 1.7 larger than that derived from ^^CO. As discussed 
above, the ^^CO is probably a better tracer of the large-scale distribution of gas, but is more 
hkely to miss material in the central region of the cloud due to uncorrected saturation. The 
results from both tracers are in reasonable agreement and we should have results from both 
in mind for comparison of the fraction of atomic gas measured here with that predicted from 
models presented in the following section. 



(2) 
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6. TIME-DEPENDENT MODELING OF THE ATOMIC TO MOLECULAR 

HYDROGEN RATIO 

6.1. Rate Equation for H2 Formation 

In order to calculate the gas phase abundances of atomic and molecular hydrogen, we 
need to model the grain surface production H2. For the present study, we assume that all 
dust grains arc spherical particles, differing only in their size. This neglects the changes in 
grain density and structure that could be a consequence of processes which initially formed 
the dust grains, or modified them by accretion or coagulation. But since the effects of grain 
evolution, particularly on the surface properties of grains which are critical for H2 formation 
and molecular depletion are so uncertain, this simplification seems a reasonable starting 
point. 

The basic equation giving the formation rate of H2 molecules on the surfaces of grains 
of a single type, e.g. grains having a specified cross section and number density (HoUenbach, 
Werner, & Salpeter 1971), is 

— -^ShI^H2''^HI < VhI > CTgrUgr . (3) 

has units of cm~^ s~^, Shi is the probability of a hydrogen atom which hits a grain 
sticking to it, and ej^j is the probability of a hydrogen molecule being formed by recombina- 
tion of two hydrogen atoms and then being desorbed from the grain^. The atomic hydrogen 
density in the gas is uhi (cm~^), the mean velocity of these atoms is < vhi > (cm s~^), 
while the number density of grains having cross section agr (cm^) is Ugr (cm^'^). 

The sticking coefficient, estimated to be 0.3 by Hollenbach & Salpeter (1971), has been 
studied in a quantum mechanical treatment by Leitch-Devlin & Williams (1985). These last 
authors considered physisorption and chemisorption on plausible grain surfaces of graphite, 
silicate, and graphite covered with a monolayer of water. The values of Shi found range 
from 0.9 to very small values. Given the uncertainties in grain composition and binding 
mechanisms, it appears reasonable to adopt the value Shi = 0.3. 

For grain surfaces with binding sites of a single type, there is an inherent competition 
between thermal desorption and H2 formation. If the binding is very strong, atoms on the 



^Shi and e/fj were multiplied together to form the H2 recombination coefficient 7 by Hollenbach, Werner, 
& Salpeter (1971), but following more recent usage we separate the two terms here. There remains some 
inconsistency in the terminology used, as e.g. Cazaux & Ticlcns (2002) use the expression "recombination 
efficiency" for 6/^2, which clearly does not include the sticking probability Shi- 



-27- 



surface will not desorb, but they will not be mobile, so the H2 formation rate will be low. On 
the other hand, if the binding is very weak, the mobility will be high, but so will be the rate 
of thermal desorption, and again, the rate of formation of H2 will be reduced. This has led 
to the concern that the formation rate of molecular hydrogen will be large only over a very 
narrow range of grain temperature, a result which would certainly be problematic given the 
wide range of temperatures characterizing clouds in which significant fractions of molecular 
hydrogen are found. 

Cazaux & Tielens (2002) circumvented this problem by postulating two types of binding 
sites on individual grains - physisorption sites with relatively weak binding, and chcmisorp- 
tion sites characterized by 10 to 100 times stronger binding. At the temperatures of dense 
molecular clouds, there will always be atoms in chemisorption sites, because the thermal 
desorption time scale far exceeds any other. Thus, any hydrogen atom hitting the grain and 
sticking will end up in a physisorption site, from which it can tunnel to a chemisorption site 
where it finds a partner. The barrier between the two types of sites is moderate (E/k = 200 
K), so that the effective mobility is high^. Cazaux & Tielens (2002) show that for grain 
temperatures between 6 K and 30 K, is unity, and drops gradually to 0.2 for a grain 
temperature of 100 K. We adopt here en^ = 1-0. 

The grain mass per unit volume implied by the single type of grain assumed in equation 
3 is given by 

where nigr is the mass of a grain, Ugas (— + nni + nue) is the total volume density of 
gas particles, < nigas > is the average mass of a particle in the gas, and GDR is the gas to 
dust ratio by mass. The average mass per gas particle is 

< f^gas >— ) \p) 

"^gas 

where nip is the mass of the proton. With this, we can write 

'^gas ^ /^\ 

UgrCTgr - ^^^^^^^ " 

The mass of a grain can be written nigr = PgrVgr, where pgr is the density of a grain and 
Vgr is its volume. Since we have no clear evidence that the grain density varies systematically 
with grain size, we adopt a single value of the grain density (typically ~ 2 gm cm~^) for grains 



^This is consistent with the most comprehensive treatment, which suggests that hydrogen atoms will be 
highly mobile under any reasonable conditions (Leitch-Devlin & Williams 1984). 
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of all sizes. Associating pgj. with the other grain parameters assumed to be independent of 
grain size, we can rewrite the H2 formation rate as 

2pgr GDR IV ■ ^ ' 

For grains of a single size, the H2 formation rate is proportional to the ratio of the grain 
cross section to its volume. We can write the H2 formation rate as 

= kH2'ngas'nHi , (8) 
where all of the other factors in equation 7 have been subsumed into 

7 _ ShI^H2 < '^HI >< nigas > CTgr . . 

~ 2pgrGDR \r ■ 

For a single "standard" grain species of radius a^, o-gr/Vgr = 3/40^. 

In Appendix A, we evaluate the effect of a grain size distribution on the rate coefficient 
for H2 formation. With the simplest assumptions that all grains have the same density, and 
that the sticking coefficient and formation efficiency are independent of grain size, the effect of 
a power law distribution of grain radii can be analyzed straightforwardly, and for the Mathis, 
Rumpl, & Nordsieck (1977; hereafter MRN) size distribution, we find the particularly simple 
result given by equation A 12 that kH2 depends only on the maximum and minimum values 
of the grain size distribution. 

Defining the total proton density in atomic and molecular hydrogen as 

no = nni + 2nH2 > (10) 

and including a standard helium abundance of 0.098 relative to hydrogen in all forms, we 

can write . 

ngas < nigas >= (1 H jnipHo = 1.39mpno , (11) 

no 



which gives 
with 



Rh2 — ^'h2 MRNl^Hino , (12) 

,/ _ SHieH2 < Vhi > 1.39mp 3 Qg .^„x 

'^H2 MRN — o7r~rrm 1^ / n n ' ' ^ ^ 

The final fraction represents the correction due to the grain size distribution, being equal 
to unity for any single grain size. Substituting standard values for the relevant parameters. 
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and normalizing to them, we find ^ 



k' = 3.5 X 10-^ ^ . "-' rnff ""^ , ' cmh-' . (14) 

y 2gcm-'i ) y Im ) \ 1.7 xW-^ cm) V "'rnaxU^mm 



Based on the discussion in Appendix A, and the values for a Mathis, Rumpl, &: Nordsieck 
(1977) grain size distribution with a^j„ = 25 A, a^ax — 10,000 A, and = 1,700 A, we find 
that the effect of the grain size distribution is to increase the formation rate coefficient by a 
factor of 3.4. We thus take a nominal value for k'jj^ — 1.2x10"^^ cm^s~^. 

The numerous factors which enter into the H2 formation rate coefficient inevitably result 
in a relatively large uncertainty in this parameter. The range of values for the sticking coef- 
ficient Shi and the H2 formation efficiency ejjj for "traditional" grains have been discussed 
earlier in this section. The grain density for grains of standard composition can vary by a 
factor ^1.5, and a similar variation in the gas to dust ratio is entirely plausible, although 
data in dark clouds is very sparse. The upper and lower limits on the grain size distribution 
enter as {amaxCi-min)~^'^ , but these quantities could easily differ from our adopted values by 
a factor ~ 2. It is also not clear whether very small grains would be effective for producing 
molecular hydrogen as their surface properties are likely very different from grains studied 
theoretically or measured in the laboratory. Combining the uncertainties in the various fac- 
tors, k'ff^ may differ from its nominal value by a factor of 5. This directly impacts the H2 
formation time scale and the steady-state HI density, both of which vary as k'jj ~^ . 



''Following almost all treatments of H2 formation, we include only the thermal velocity of the hydrogen 
atoms in equation In principle, wc might also consider the nonthermal component of the line widths, which 
is generally considered to be turbulent in nature. However, this will affect processes such as the coUsion of a 
hydrogen atom with a grain only if the characteristic size scale of the turbulence is comparable to or smaller 
than the mean free path of the atoms. Taking the most conservative approach including only collisions 
between atoms and grains, this is approximately lO^^/ugas cm. This very small distance is far smaller than 
the scale size on which line widths in cloud cores approach their purely thermal values (Barranco & Goodman 
1998; Goodman et al. 1998). Thus it appears justified to include only the thermal velocity in equation 
14., which gives the H2 formation rate varying as T'^ "''. Fortunately, the nonthermal contribution to the 
linewidth for HINSA is, unlike the situation for HI emission and most absorption, quite modest compared 
to the thermal linewidth (as seen in Table 2 of Li & Goldsmith (2003) and Table 6 here) and thus the 
uncertainty introduced by considering only the thermal velocity for the formation rate is less than a factor 
of 2. 



-30- 



6.2. Time Dependence of HI and H2 Densities 



An idealized scenario for the evolution of the abundance of HI and H2 consists of a 
cloud in which the hydrogen is initially entirely in atomic form, maintained in this state by 
photodestruction. At t = the extinction is suddenly increased so that photodestruction in 
the bulk of the cloud volume (regions with visual extinction greater than ~ 0.5 mag) can 
be neglected. This is justified since the photodissociation rate of H2 drops dramatically as 
a function of increasing column from the cloud surface, due to effective self-shielding (e.g. 
van Dishoeck & Black 1988; Draine & Bertoldi 1996). 

To model the evolution of this cloud, we assume that the only destruction pathway for 
molecular hydrogen is cosmic ray ionization ^ , which is assumed to occur at a rate (^Hq s"^. 
The value of in dense clouds has been found to he in the range 10~^^ to 10~^^ s~^ (CaseUi, 
Walmsley, Terzieva, & Herbst 1998), but in other studies has been more narrowly defined 
to be equal to 5 x IQ-^'^ s'^ (Bergin et al. 1999), and (5.2 ± 3.6) x 10"^^ s"^ (van der Tak 
& van Diskhoeck 2000). This is consistent with results of Doty et al. (2002) from detailed 
modeling of the source AFGL2591. 

A much higher ionization rate, = 1-2 x 10~^^ s~^ found in diffuse clouds (McCall 
et al. 2003), can possibly be reconciled with the above if there is a component of the cosmic 
ray spectrum which penetrates diffuse but not dense clouds. Doty, Schoier, & van Dishoeck 
(2004) suggest that the much higher rate which they derive for IRAS 16293-2422 is actually 
being dominated by X-ray ionization from the central star, and that the true cosmic ray 
ionization rate is much lower. We here adopt the average of the values obtained by van der 
Tak & van Dishoeck, (h^ = 5.2x10"^^ s-^ 

We assume that the evolution of the cloud takes place at a constant density, and write the 
equation expressing the time dependence of the molecular hydrogen (denoting the formation 
rate coefficient simply by k') as 



Defining the fractional abundance Xi (of atomic and molecular hydrogen) as the density of 
species i divided by the total proton density, we find that 



dnH2 
dt 



(15) 




— = k xiuo - CH2X2 ■ 



(16) 



^There are cases in the literature where it is unclear whether one is dealing with the cosmic ray ionization 
rate per H atom or the rate per H2 molecule (in molecular clouds). Since the regions studied here are largely 
molecular, we adopt the latter definition. 
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Substituting equation 10 we can rewrite this as 

^ = k'no - {2k'no + Ch,)x2 ■ (17) 

The solution is the time dependence of the fractional abundance of molecular hydrogen, 
given by 

= [1 - exp{-t/THi^H,)] ■ (18) 

2k' no + 

The fractional abundance of atomic hydrogen is 1 - 2x2, and hence is given by 

x.it) = 1 - J^'""' [1 - expi-t/rni^H,)] ■ (19) 
The time constant for HI to H2 conversion is given by 

The steady-state molecular and atomic gas fractions can be found from equations 18 and 
19, but in the limit of present interest k'no 3> Ch2 we obtain 

= 2k ■ '''' 

The steady-state values for the constituents of the gas are then 
and 

(h2 Ch2 /r,QN 

xi —— or uhi^^ . (23) 
2k'no 2k' 

This limit thus corresponds to our having essentially molecular gas, but the implication 
of equation 23 is that there is a constant density of atomic gas given by 

UHiit » TH1-.H2) =^Hi^% ■ (24) 

In Figure 10 we show the evolution of the atomic hydrogen fractional abundance as a 
function of time, for clouds in which the total proton density is fixed. From the discussion in 
Section 6.1 and Appendix A, we have adopted k' — 1.2 x 10"^''' cm^s~^ and C/Za = ^-^ x 10"^^ 
s~^. This yields an equilibrium atomic hydrogen density n*jjj equal to 2.2 cm~^. Equation 
24 is not at all new; it was given by Solomon & Werner (1971), but their use of a much 
higher value of C,h2 resulted in an estimate of the equilibrium abundance of atomic hydrogen 
equal to 50 cm~^, a value clearly ruled out by our data. 
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Fig. 10. — Time dependence of the fractional abundance of atomic hydrogen for clouds 
of fixed total proton density (in cm~'^). The curves of different total proton densities are 
logarithmically spaced at intervals of 0.25. The results for five sources studied here, located in 
terms of the atomic hydrogen to total proton density ratio and total proton density, are also 
included. The four sources with satisfactory solutions are indicated by open circles. B227, 
for which the fractional HI abundance is very close to the steady state value, is indicated by 
an arrowed box at the youngest age reasonably consistent with the data. 
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7. DISCUSSION AND IMPLICATIONS FOR CLOUD EVOLUTION 

7.1. Cloud Age Based on Time— Dependent HI Fraction 

For the nominal values we have adopted, the simple form of thi^H2 given by equation 
21 will apply for uq » n* = 2.2 cm~^. In this hmit, which is certainly relevant for dense 
clouds, the time constant is 

8.3 X 10^6 2.6 X 10^ 

THI-.H2 = s = yr . (25) 

no no 

Assuming that the clouds we have studied are characterized by densities determined from 
the ^^CO observations, no = 2 to 6 xlO^ cm~^, and thi^H2 = 0-5 to 1.3 xlO^ yr. Using the 
larger densities from C-'^^O reduces this characteristic timescale by a factor ~ two ^ . As seen 
in Figure 10, the fact that the clouds observed have HI fractional abundances relatively close 
to steady-state values immediately indicates that the time elapsed since initiation of HI to H2 
conversion must be several times greater than thi^H2 ot at least a few million yr. This "age" 
applies to the bulk of the molecular cloud as traced in ^^CO and HINSA ^'^ . For four of the 
five clouds studied, we obtain a unique solution for the elapsed time. This quantity, which 
we denote Tdoud, is defined by the total proton density (twice the molecular hydrogen density 
obtained from the ^'^CO column density and source size), and the fractional abundance of 
atomic hydrogen. While keeping in mind the uncertainty of the "age" as determined here, 
for these four dark clouds, this quantity lies in the range 10^-^ yr < Tdoud < 10^ yr. The 
fractional abundance of atomic hydrogen in B227 is almost exactly the steady state value, 
and choosing a minimum time when this quantity will have dropped to close to this value 
leads to Tdoud > 10''' yr. 

These relatively great cloud ages would pose a problem if there were clear evidence 
that these clouds were dynamically evolving. To evaluate this possibility, we have used the 
gaussian density distribution and the parameters derived above. Using the results from 
Appendix B and the masses and line widths given in Table 6, we calculate the potential and 
kinetic energies, ignoring surface effects, the external medium, and any contribution of the 
magnetic field. The clear result is that these clouds are close to virial equilibrium, with the 



^The central densities derived from our observations do not reflect the true central density of the cloud 
because the C^^O J = 1 ^ transition is not sensitive to densities > 10^ cm~^, and we have smoothed the 
FCRAO data to 3'. 3 beamwidth for comparison with the Arecibo HI data. 

Again, this does not apply to the cloud cores, particularly that of L1544 where densities are 10^ - 10® 
cm~^ (Caselli et al. 1999; Tafalla et al. 2002) and where the HI — > H2 conversion timescale will be much 
shorter. 
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possible exception of L1544, for which the numbers are particularly uncertain due to the 
incompletely defined extent of the cloud, but which appears to have —U noticeably larger 
than 2T. While this does not define the evolutionary path these clouds have followed, it 
does suggest that it is entirely plausible that they have spent several million years evolving 
to their present state^^ . 

In fact, we feel that these vahics for cloud ages must, for a number of reasons, be 
considered as lower limits. First, it is unlikely that the evolution of these clouds proceeds 
at constant density. Typical models of interstellar cloud evolution indicate that the path 
followed starts with lower-density atomic clouds and proceeds to higher density. This is not a 
requirement, in that if, for example, a compressive event initiates the process by compressing 
the gas and increasing the visual extinction, the density could immediately be increased to 
something like its final value, but in fact, no atomic clouds with even the moderate densities 
given in Table 5 have been observed. If cloud evolution includes an increase in the density 
occurring simultaneously with the HI to H2 conversion, then much of the evolution will 
have taken place at lower densities than those now observed, and the time scale will be 
correspondingly longer. However, if the clouds were much more turbulent in the past and 
the turbulence were on a scale that affected the velocities of hydrogen atoms colliding with 
grains (see footnote 7), then the lower collision rate and H2 formation rate resulting from 
the lower density would be partially offset by the higher Hl-grain velocities. 

Second, we may consider what happens if the column density is not sufficient to make 
photodestruction of H2 insignificant in the central regions of these clouds. To first order, we 
can consider photodestruction of H2 simply as adding to the cosmic ray destruction rate in 
equation 15. While formally increasing the time constant (equation 20), for clouds of high 
density, thi-^H2 is unaffected, as long as uq » [Ccr + CpHOTol/'^k' . The models of Draine & 
Bertoldi (1996) suggest that the H2 photodestruction rate will be comparable to the cosmic 
ray destruction rate adopted here for a column density of H2 ~ 10^^ cm~^. Thus, while 
there could be a modest contribution from photodissociation to the total H2 destruction 
rate, it will not be enough to change the time constant. However, the photodestruction will 
increase the steady-state abundance of HI. The fact that the values we have measured are 
only modestly above those expected solely from cosmic ray destruction of H2 indicates that 
the effect is modest. If photodestruction plays a significant role, our solutions become lower 
limits to Tcioud since they would, in effect, be steady-state values. 



^^Caselli et al. (1999) suggest that the core of L1544 is collapsing with a much shorter characteristic time 
scale. Only a very small fraction of the cloud is involved in the collapse, suggesting that the evolution of 
this region is, in fact, decoupled from that of the bulk of the cloud which we are studying here. 
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Third, as discussed in Section 5, we may have underestimated the column densities of 
molecular hydrogen, and hence ng, along the lines of sight to the cores of these clouds. This 
may result in our having modestly overestimated whi/wq in these directions. This would 
have the effect of making the clouds look "younger" than they actually are. 

Fourth, we must consider the possibility that there is not a reservoir of chemisorbed 
hydrogen atoms on the grain as postulated by (Cazaux & Tielens 2002, 2004). Then, if 
the gas phase HI density and thus the adsorption rate is sufficiently low, one can reach a 
situation in which <^hi>^ the average number of hydrogen atoms on the grain, is less than 
unity. A newly-adsorbed incident atom may not find another atom with which to combine, 
and the H2 formation rate will be reduced. This issue has been discussed by Biham et al. 
(1998) and Katz et al. (1999), and is consistent with laboratory measurements (e.g. Vidali et 
al. 1998). There is some concern whether the laboratory measurements are really applicable 
to the surfaces of interstellar grains, so that it seems appropriate to be cautious about the 
relevance of the measurements. If this picture is correct, the rate of H2 formation for small 
grains is also reduced since <Hhi> will also be less than unity. These effects combine 
to reduce the overall H2 formation rate, particularly at relatively late times in the cloud 
evolution, when uhi < 100 cm~^. 

Additional support for our deteriming the minimum age of dark molecular clouds comes 
from modeling of the thermal evolution of gas which is in the process of converting hydrogen 
from atomic to molecular form. Each hydrogen molecule formed releases a few cV of energy 
into the gas. The exact value is uncertain due to unknown distribution of the 4.5 cV molecular 
binding energy between grain heating, internal energy of the H2 (which is radiated from the 
cloud), and kinetic energy of the H2, which eventually is thermalized (Flower and Pineau des 
Forets 1990). Under steady state conditions and high densities as we have seen in Section 
6.2, the density of atomic hydrogen is ~ few cm~^. This density of HI results from the 
balance between formation of HI via cosmic ray destruction of H2 and destruction of HI by 
formation of HI on grains. In this situation, the heating from H2 formation can be considered 
as a part of the cosmic ray heating, as it is the latter which supplies the energy to break the 
molecular bonds. 

Early in the evolution of a cloud from atomic to molecular form, if we assume that 
the beginning of the process is characterized by a sharp increase in density to something 
approximating its current value, then the heating rate from H2 formation is much larger, by 
a factor no/n|^j ~ 10'^ (Goldsmith & Langer 1978). The conversion of atomic to molecular 
hydrogen thus results in a reheating of the cloud (which had been cooling after the initial 
compression), to a temperature ~ 100 K (Flower and Pineau des Forets 1990). This re- 
heating phase lasts for a time approximately equal to the characteristic HI H2 conversion 
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time (equation 20), which is plausibly on the order of a million years. The fact that we 
see dark clouds which have obviously cooled to temperatures an order of magnitude below 
that sustained by the initial formation of molecular hydrogen is a further indication that the 
initiation of the molecular cloud phase must have taken place a minimum of a few million 
years in the past. 

7.2. Other Effects and Concerns 

Although the time-dependent evolution of the HI fractional abundance in the scenario 
developed above seems plausible, there are a number of other processes which may affect 
the abundance of atomic hydrogen and thus confuse this picture. One such effect is turbu- 
lent diffusion, in which the enhanced ability of material to move over significant distances 
increases the effect of concentration gradients (Xie, Allen, & Langer 1995). Turbulent dif- 
fusion has been apphed by Willacy, Langer, & Allen (2002) to the case of atomic hydrogen 
in dark clouds. Their results suggest that a quite modest values of the turbulent diffusion 
constant, K < 5 x 10^^ cm^ s~^ are sufficient to produce atomic hydrogen densities in the 
centers of clouds similar to what we measure. This is not surprising, in the sense that our 
values are equal to or modestly greater than those expected for steady-state conditions with 
no turbulent diffusion. 

It appears that while these values of K may be physically plausible, they are significantly 
smaller than those required for this model to reproduced various other chemical abundances 
in dark clouds. It is also the case that the atomic hydrogen column densities predicted by 
this model are significantly in excess of those that we observe. Thus, it is premature to say 
that the present observations indicate, in any significant way, that turbulent diffusion plays 
a very significant role. The results of Willacy, Langer, & Allen (2002) for the density of HI 
in the center of the model cloud are for 5 x 10^ yr of evolution of the cloud, and it is plausible 
that the timescale we have derived above still is relevant for reducing the HI density in the 
centers of dark clouds to the low values we have obtained. 

A related model suggests large-scale circulation of gas in the cloud (Chieze & Pineau 
des Forets 1989). In this model, parcels of gas are exchanged between the outer, relatively 
unshielded portions of the cloud, and the interior. These authors analyzed the effect on a 
number of molecular species, but did not consider HI. However, it is reasonable that the 
effect would be to increase the HI density inside the cloud, as is the case for other species 
which are typically abundant only in regions of small visual extinction. While suggestive, it 
is difficult to make any quantitative comparisons with this model as no results for atomic 
hydrogen are reported. 
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We should consider the correlation between cold HI and molecular column density pre- 
sented earlier in Figure 9. It is evident that in the context of a steady-state for HI production 
and destruction, the HI column density is proportional to the line of sight extent of the cloud, 
since the number density of atomic hydrogen is a constant. In this sense it is very different 
than the column density of e.g. ^^CO, which is the integral of the number density of that 
species along the line of sight. For the HI, the extent of the cloud is be defined by two 
factors. The first of these is the size of the region in which the HI number density is more 
or less constant. Evidently, the fractional abundance of the HI will increase dramatically in 
the portion of the cloud with low extinction, but outside of some radius, the number density 
will become insignificant. The second factor is basically observational: since the absorption 
coefficient for HI varies inversely as the temperature, atomic hydrogen in the outer layers of 
the cloud is much less detectable than that in the cold interior. This effect is significant as 
the outer layers of the cloud even with standard radiation field can be at temperatures ~ 
100 K, an order of magnitude higher than those in the cloud center. 

It is appropriate here to follow up on this issue of warm HI surrounding that which we 
feel is well-mixed with the molecular material. Models of basically static clouds embedded 
in a much warmer, diffuse interstellar medium and irradiated by an interstellar radiation 
field (ISRF) have envelopes ( "onion skins" ) of ionized hydrogen, and neutral hydrogen, and 
ionized carbon, and neutral carbon (e.g. Abgrall et al. 1992; Le Bourlot et al. 1993). The 
self-shielding of the H2 results in its abundance rising when the visible optical depth Ty 
reaches 10~^, and becoming comparable to that of HI for = 10~^. The temperature will 
be 30 K and 100 K, for enhancement factors of the radiation field relative to the standard 
ISRF of 1 and 10, respectively. 

The HI halos seen in emission around molecular clouds (e.g. Wannier, Lichten, & Morris 
1983; Andersson, Roger, & Wannier 1992) correspond to these warm cloud edges where the 
fractional abundance of H2 is dropping significantly. However, as discussed earlier, the very 
narrow line widths and low minimum line intensities that we are dealing with here indicate 
that we are observing much cooler and also less turbulent gas. The smaller spatial extent of 
HINSA absorption compared to the CI emission as well as its narrower line width (Section 
3.7) is additional strong evidence that the HINSA column density is not significantly affected 
by the warm neutral/ionized cloud halos . 

The cold HI and molecular column densities, as indicated by Figure 9 (as well as by the 
maps of their distributions), are quite well correlated. Given the very different distributions 



^^This is the same conclusion reached by Flynn et al. (2002) based on modeling the HI self absorption 
produced by cold clouds in the Galactic disk. 
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expected from theoretical models, this appears somewhat surprising, but may reflect to some 
degree the basic geometry of the clouds. While there is clearly a significant density gradient, 
these regions have a basically convex geometry in three dimensions, so that lines of sight with 
large column density do encompass larger dimensions with sufficient material to provide a 
significant contribution to the column of cold HI. More complete models including the spatial 
and temporal evolution of atomic and molecular hydrogen, as well as other species, will be 
necessary to assess how accurately the atomic to molecular abundance ratio in the centers 
of these clouds can be used to determine their age. 



7.3. Upper Limit on the Cosmic Ray Ionization Rate 

We can utilize our measurement of the atomic hydrogen density to set a useful upper 
limit on the cosmic ray ionization rate in reasonably well-shielded dense clouds. From 
equation 24 we can write the cosmic ray ionization rate in terms of the steady state atomic 
hydrogen density and the H2 formation rate as 

Ch, = 2k'H,n*Hi . (26) 

The HI density we have measured is an upper limit to the steady state density of atomic 
hydrogen, since processes such as turbulent diffusion, imperfect shielding against photodis- 
sociation, and time-dependent cloud evolution all work to increase the atomic hydrogen 
density. As discussed in Section 6.1, the uncertainties in k^^ ^-^^ substantial, but an upper 
hmit to k^2 equal to 6xl0~^^ cm^ s"^ is reasonable. If we interpret our observations as 
giving uhi — 4 cm~^ as the upper hmit to n*jjj, we find that C,H2 < 5 x 10~^^ s~^. 

Our data, in common with many other observations of dark clouds, constrain the cosmic 
ray fiux through the heating provided to dark clouds. In regions with visual extinction greater 
than a few magitudes, cosmic rays are expected to be the dominant heating source for the 
interstellar gas. At densities < 10^ cm~^ gas-dust couphng has only a very minor effect on 
the gas temperature, which is determined by the balance between cosmic ray heating and 
molecular line cooling (Goldsmith 2001). The latter varies as the gas temperature to the 
2.4 power for n^jj = 10^ cm~^ and 2.7 power for n^j = 10^ cm~^ 

The minimum gas temperature is obtained with no depletion of coolant species, and for 
n^2 in the range 10^ to 10^ cm~^, we find that Tgas — Tqq^-^, where g is the scaling factor 
for the cosmic ray heating rate relative to a reference value 10~^^ erg cm^^ s~^. Tq is 
equal to 10 K for n^f^ = 10^ cm~^ and 13 K for n^^ = 10^ cm~^. 

For our adopted cosmic ray ionization rate C(H2) = 5.2x10"^^ s~^, and taking AQ = 7 
eV as the average energy transferred to the gas per cosmic ray ionization in a region of low 
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fractional ionization (Cravens & Dalgano 1978), we obtain a heating rate F = 5.8xl0~^^nij2 
erg cm~^ s~^. This heating rate corresponds to g = 0.58, which yields gas temperatures Tgas 
= 8.0 K for n//2 = 10^ cm~^ and 10.5 K for 11^2 = 10"^ cm^"^. These are in good agreement 
with the temperatures of clouds, 8 K to 10 K determined in many studies using ^^CO (e.g. 
Li & Goldsmith 2003), and ammonia (Tafalla et al. 2002). These temperatures are also in 
excellent agreement with the low central temperatures inferred from modeling preprotostellar 
cores (Evans et al. 2001). 

To be consistent with measured values of the kinetic temperature, wc cannot have g 
substantially greater than unity. Taking a reasonable upper limit to be g = 1 yields Tcr < 
lxlO~^''nj:/2 6rg cm~^ s~^, and gives gas temperatures Tgas = 10 K for nj^j = 10^ cm~^ and 
13 K for = 10^ cm~^. A cosmic ray heating rate significantly greater than this value 
thus seems unlikely to characterize the dense regions in dark clouds. The much higher values 
that are sometimes inferred (e.g. McCall et al. 2003; Liszt 2003) may be the result of local 
conditions or of more diffuse compared to denser regions, since it is relatively more difficult 
to depress the heating rate than to enhance it. The low temperature of molecular clouds 
such as those observed here also suggests that the upper limit for the cosmic ray ionization 
rate derived above may well be overly generous and that the actual value of does not 
significantly exceed 10^^^ s"-*-. 



8. SUMMARY 

We have carried out a study of four dark interstellar clouds to define better the location 
of the cold atomic hydrogen relative to the molecular component within these regions. The 
HI is traced through observations of narrow self absorption features (HINSA), which are seen 
against the general Galactic background. The molecular clouds are traced by observations 
of OH (Arecibo) and ^^CO and C^^O obtained at FCRAO. For comparisons between atomic 
and molecular species, our results are limited to regions > 3' in size, and thus refer to 
moderate-scale structure of these clouds. The most quiescent well-shielded central regions, 
in which molecular depletion may be significant, are heavily beam-diluted in this study and 
thus are unlikely to be affecting our conclusions. 

The extent of the cold HI is close to, but somewhat smaller than that of the ^^CO, and 
larger than the size of the C^^O emitting region. This is consistent with results previously 
obtained by Li & Goldsmith (2003), in which the nonthermal line widths of the species 
follow this same ordering. The cohimn density of cold HI is moderately well correlated with 
that of ^"^CO, which is consistent with the general impressions of the maps of the different 
species, and which further indicates that the atomic hydrogen we are observing is not in outer 
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envelopes of these clouds, but is predominantly in well-shielded cold regions, again consistent 
with the low temperatures for the HINSA-producing gas derived by Li & Goldsmith (2003). 

The mapping of ^^CO and C^^O isotopologues reveals that the ratio of the column 
density of these species rises sharply at the edges of the clouds. This is unlikely to be a result 
of any radiative transfer effect, but rather suggests that the abundance ratio is enhanced in 
regions of low extinction. This effect is consistent with chemical isotopic fractionation, which 
has been predicted to be significant where the abundance of ionized carbon is significant, 
but the temperature is still quite low. The enhancement of the ^^C in carbon monoxide 
offsets the drop in the abundance of this molecule in the outer regions of clouds due to 
photodestruction. The result is that the size of the region over which the abundance of ^'^CO 
relative to H2 remains relatively uniform is increased, and we feel that this isotopologue is 
to be preferred as a tracer of the large-scale cloud structure, although saturation effects in 
regions of largest column density must be considered. 

Using the HINSA and ^^CO data, we model these clouds and obtain the central densities 
of atomic hydrogen and H2. The central densities of H2 are 800 to ~ 3000 cm~^, and those 
of cold HI are between 2 and 6 cm~^. The ratio of the atomic hydrogen to total proton 
densities ranges from 0.5 to 3 xlO~^. The values of thh^ from C^^O are a factor of two larger 
than those from ^^CO, and the fractional abundances of HI a corresponding factor lower. 
These absolute and fractional HI abundances are close to, but somewhat larger than those 
expected for steady-state balance between HI production by cosmic ray ionization of H2, 
and production of this molecule by grain surface recombination of HI. The effect of a power 
law grain size distribution on the formation rate of H2 molecules has been evaluated, and for 
a MRN distribution with reasonable upper and lower limits, we find that that H2 production 
rate is increased by a factor of 3.4 over that for standard grains of radius 1700 A. 

A number of processes including mass exchange and turbulent diffusion may be playing 
a role in determining the abundance of atomic hydrogen in dark clouds, but it is not certain 
whether their contribution is significant. We assume that these clouds have been evolving at 
constant density since the time that the atomic to molecular transition was initiated by some 
process which suddently increased the visual extinction and eliminating photodestruction 
of H2. In this idealized picture, the HI to H2 ratio drops monotonically with time, and 
the time constant for the conversion is given by thi-^H2 = 2.7 x 10^/no yr, where no is 
the total proton density. For the central densities of these sources, the time constants 
are ~ 10*^ yr, and the ages of these clouds are 3 to 10 million years. The relatively low 
fractional abundance of atomic hydrogen makes it extremely unlikely that the time since the 
compressive event starting the atomic to molecular transition is significantly less than this, 
since photodestruction and the other processes mentioned above increase the HI abundance. 
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and thereby require that the grain surface production of H2 be quite close to completion. 
The conversion of atomic to molecular hydrogen sets a quite large lower limit to the time 
scale for the "molecularization" of dense clouds, and thus for the overall process of star 
formation. 
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A. Effect of Grain Size Distribution on H2 Formation Rate 

In order to investigate the effect of a grain size distribution, we start from equation 
3. For a grain size distribution, again assuming other quantities are independent of grain 
size, we use equation 9, but must take the appropriate grain cross section and grain volume 
averaged over the grain size distribution . We can write the formation rate coefficient for 
a grain size distribution as 

" 2p,rGDR V,r ■ ^ ^ 

Assuming a distribution of spherical grains, and taking the relative number of grains 
having radius between a and a + da to be n{a)da, we have 



J n{a)a{a)da , (A2) 



and 

Vgr = I n{a)V{a)da . (A3) 



^"^Thc quantity which needs to exphcitly be integrated over a is the grain cross section. Integrating the 
grain volume over a is necessary only to normalize by a fixed total mass in grains, given that we have assumed 
that all grains have the same density. 
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We assume a power law grain size distribution having the form 

n{a) = n{amin){-^)~'' , (A4) 

where a^i„ is the minimum grain radius and Umax is the maximum grain radius (a^j„ < a < 
dmax)- The power law exponent, a, is 3.5 in MRN model, but is here left as a free parameter. 
The total density of grains (cm~^) is given by^^ 

n{a)da= / n{amin){ y^da . (A5) 

Evaluating equation A5, we find 

n{amin) = _ -, , (A6) 

where R is the ratio of maximum to minimum grain radius: 

. (A7) 

The ratio of the integrated grain cross section to integrated grain volume is 

"4 - a 



Sgj. 3 

Vgr ^(^miri 



3 — a 



The effect of grain size distributions characterized by different values of a on T,gr/Vgr is 
shown as a function of R in Figure 11. 

In the limit of no grain size distribution (R — > 1) we recover the expected result for 
spherical grains of a single radius, a^. 



y 



We can consider various values of amax and amin by treating amin and R as free parameters, 
but it is helpful to normalize to the cross section to volume ratio for a single grain size, which 
gives us 

- 1 (AlO) 



Sgj, / 3 
Vgr I 4a5 



Equations A8 and AlO are general, but of particular relevance to suggested astronomical 
grain size distributions, these expressions are well-behaved for all values of a > 0, including 
e.g. 3 and 4. 



^^The dimensions of n(a) and n{amin) are evidently cm 
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log(R) = log(a„,ya„,.J 



Fig. 11. — Effect of grain size distribution with n{a) oc or" on tfie ratio of tlie integrated grain 
cross section to tlie integrated grain volume. The H2 formation rate is directly proportional 
to this ratio. The vertical axis is normalized to the value of E^^/V^r for grains of a single 
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For a fixed value of Umin, ^gr/Vgr decreases as R increases. For i? » 1, the ratio varies 
as for < a < 2.5, and has a logarithmic dependence for R = 3.0 and 4.0. The MRN 
grain size distribution is characterized by a = 3.5. In this case, we find that, for all R, 



Tigr / 3 

Vgr I 4:Clfnin 



= . (All) 

o=3.5 



This particularly simple result facilitates assessing the impact of this grain size distribution. 
The unnormalized ratio of the integrated grain cross section to integrated grain volume can 
be written 

^ = ^ ^ , (A12) 

^gr a=3.5 ^ '\/^max ^min 

indicating that for this value of a, the H2 formation rate is inversely proportional to the 
geometric mean of the extremes of the grain size distribution. 

In order to interpret the effect of the grain size distribution quantitatively, we need to 
consider specific values of the "standard" single grain size, a^, along with the specifics of 
the GSD defined by amin and R, in addition to a. Normalizing relative to the single-sized 
grains, we find 

^ /— = = . (A13) 

The ratio of integrated grain cross section to integrated grain volume relative to the ratio 
for grains of a single radius is thus determined by the ratio of that single grain radius to 
the geometric mean of the maximum and minimum grain radii. A smaller grain has a larger 
surface to volume ratio, so a grain size distribution including grains that are all smaller 
than the standard grain, will have E^^/V^r larger than 3/4as. It is more meaningful to 
compare grain size distributions with the constraint amin < Os < Cimax, so that for i? — > 1, 
{12gr/Vgr)/{3/4as) O's/O'min — 1- Notc that as long as amin < < amax, the normalized 
grain cross section per unit grain volume is independent of the reference grain size a^. 

There is a vast and ever-expanding literature on the subject of the range of grain sizes 
in interstellar clouds. These should be considered relative to the "standard" dust grain size, 
which is based on fitting behavior of the extinction in the visible region of the spectrum 
(Spitzer 1978, ch 7). A representative value of is 1700 A. 

The more diverse dust grain population includes grains having radii as small as 10 
A which are transiently heated to ^ 10^ K in refiection nebulae by absorption of a single 
photon (Sellgren 1984). The infrared emission of small grains is also seen from the edges 
of molecular clouds exposed to fairly standard interstellar radiation field (Beichman et al. 
1988). The need for small grains (radius 5-20 A) to explain diffuse infrared emission in 



-45 - 



a number of galaxies in addition to the Milky Way was emphasized by Rowan-Robinson 
(1992). Specific spectral features seen in the infrared led to recognition that the population 
of small grains may include polycyclic aromatic hydrocarbon (PAH) molecules (e.g. Leger & 
Puget 1984; Omont 1986; Ryter, Puget, & Perauh 1987; Puget & Lcgcr 1989), although 
very small grains may well have a significant graphite component (Desert et al. 1986). 

There has been increasing evidence that a larger-than-standard grain size component 
exists in dense regions of the interstellar medium. The earliest indications of "grain growth" 
in these regions came from visual/infrared photometry (Carrasco, Strom, & Strom 1973) 
and the wavelength dependence of polarization (Carrasco, Strom, & Strom 1973; Brown 
& Zuckerman 1975; Serkowski, Mathewson, & Ford 1975; Whittet & van Breda 1978; 
Breger, Gehrz, & Hackwell 1981). Large grains (radius ~ 5000 A) have been invoked to 
reproduce the spectral properties of reflection nebulae (Pendleton, Tielens, & Werner 1990) 
and 30 //m (3x10^ A) radius grains added to a dust model to explain the emission observed 
at miUimeter wavelengths (Rowan-Robinson 1992). 

Grain growth has also been inferred from the dependence of the long-wavelength (mil- 
limeter and submillimeter) spectral index of the dust emissivity on the column density and 
volume density in molecular clouds (Woody et al. 1989; Walker, Adams, & Lada 1990; 
Goldsmith, Bergin, & Lis 1997; Visser et al. 1998). These observations indicate that the 
frequency dependence (characterized by power law exponent 0) of the grain emissivity or 
optical depth in dense regions of molecular clouds has f3 less than 2, the standard value 
characteristic of grains in the general interstellar medium (e.g. Kogut et al. 1996). The 
observations of Stepnik et al. (2003) suggest that grain coagulation in a dense filament in 
the Taurus complex has resulted in a deficiency of small grains together with an enhance- 
ment of the submillimeter wavelength emissivity. Grain coagulation is a very complex and 
imperfectly understood process (see Draine 1985 and references therein), but can result in 
fiuffy and highly non-spherical grains (Wright 1987). The MRN grain size distribution may 
be significantly modified in regions where dsut coagulation has taken place, with a resulting 
uncertain effect on the H2 formation rate. 

Evidence for grain growth is very strong in even denser regions, particularly protostellar 
disks, which suggests that this is, in fact, a continuous process throughout all phases of the 
dense interstellar medium. In choosing nominal values for the maximum and minum grain 
sizes for use in evaluating the H2 formation rate in fairly standard clouds, there is certainly 
a large range of possible values, but the effect of the uncertainty is mitigated by the fact 
that in e.g. equation A13 the dependence is on the geometric mean of a^ax and a^j„. We 
adopt nominal values for radii amax = 10,000 A and amin = 25 A, which together with the 
the standard grain radius Og = 1,700 A result in the grain size distribution increasing the H2 
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formation rate by a factor of 3.4. This modest but significant value suggests that a plausible 
grain size distribution does increase the H2 formation rate, but only in extreme cases where 
the grain size distribution has been heavily skewed to cither very large or very small grains, 
will the formation rate coefficient k'jj^ be drastically affected by the grain size distribution. 



B. Potential Energy of Cloud with Gaussian Density Distribution 

If we take the density distribution for a spherically symmetric cloud to be given by 

n(r) = noe-('^/''°)' , (Bl) 

the mass contained within radius r is given by 

M(r) = ^^rl[^/^^ erfc{—) - 2(— )e"^^^'] , (B2) 

ro ro 

where erfc is the error function. The total mass of the cloud is then given by 

M = TT^/^no < nigas > 4 , (B3) 
where < nigas > is the average mass of a gas particle, defined in equation 5. 
The potential energy is given by the integral 

U=-G r^Wr) 
Jo r 

where dm{r) is the mass in a shell at radius r. The integral is well-defined, and gives us the 
result 

U = -'^. (B5) 



A fraction 0.90 of the mass is contained within radius rgo = 1.768ro, and expressed in terms 
of this parameter, U — — 0.71 GM^ /rgo, very similar to that for a r~^ density distribution. 
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Table 1. Source Parameters 



Name 


RA(2000) 


Decl (2000) 


Distance 
(pc) 


Reference 


LI 544 


05 04 18.1 


25 11 08 


140 


Elias (1978) 




B227 


06 07 28.4 


19 28 04 


400 


Bok & McCarthy (1974); 


Arquilla (1984) 










Tomita, Saito, & Ohtani 


(1979) 


L1574 


06 08 05.0 


18 28 12 


300 


Kawamura et al. (1998) 




CB45 


06 08 45.9 


17 53 15 


300 


Kawamura et al. (1998) 





Table 2. Peak T^,, (K) at Reference 
Positions of Sources 



Name 


OH 






L1544(0,0) ^ 


1.5 


5.9 


2.4 


B227(0,0) 


0.36 


2.6 


0.41 


L1574(0,9) ^ 


0.46 


2.3 


0.44 


CB45(-3,18) 


0.44 


3.4 


0.79 



^offsets in arcminutes 
^peak of blue emission 
•^stronger peak in the cloud 
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Table 3. Column Densities at Reference Positions 



Name 


N(HINSA) 


N(OH) 


N(i3C0) 


N(Ci«0) 


N(H2)^ 


N(H2)'^ 




10^9 cm-2 


10^^ cm-2 


10^^ cm-2 


10^^ cm-2 


10^1 cm-2 


10^1 cm-2 


L1544 


.42 


1.7 


70^. 


12. 


6.0 


8.4 


B227 


.61 


.80 


31. 


2.6 


3.2 


2.8 


L1574r (0,0) 


.83 


2.3 


9.4 




1.7 




L1574b (0,9) 


1.5 


2.1 


36. 


5.1 


3.6 


4.3 


CB45 (-3,18) 


2.7 


1.0 


62. 


6.8 


5.4 


5.3 



'^From ^'^CO (column 4) and expression on p 603 of Frerking, Langer, & Wilson (1982) for 
i^CO in Taurus 

''From C^^O (column 5) and expression on p 603 of Frerking, Langer, & Wilson (1982) for 
C^^O in Taurus 

'^Corrected for saturation as described in the text 
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Table 4. Cloud Dimensions and Central Densities of Different Species 



Name 


Species 




AzpWHM LOS^ 


^cen 











(10^** cm) 


(cm 




L1544 


HINSA 


>11.5 


>1.3 


<3.0 




T,1 1^44 


13pQ 




^ i-.O 




-3 








1 9 


Q 4v 1 D" 


-4 


B227 


HINSA 


8.1 


2.9 


2.0 




B227 


i^CO 


8.4 


3.0 


9.6x10- 


-4 


B227 


Ci«0 


5.2 


1.9 


1.3x10- 


-4 


L1574r 


HINSA 


7.4 


2.0 


3.9 




L1574r 


13C0 


7.7 


2.1 


4.2x10- 


-4 


L 1574b 


HINSA 


8.8 


2.4 


5.9 




L1574b 


i^CO 


11.2 


3.0 


1.1x10" 


-3 


L1574b 


ci«o 


6.2 


1.7 


2.8x10- 


-4 


CB45 


HINSA 


15.6 


4.2 


6.0 




CB45 


I'^CO 


12.0 


3.0 


1.9x10- 


-3 


CB45 


c^^o 


10.0 


2.7 


2.4x10- 


-4 



^Geometric mean of two dimensions in the plane of the sky 
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Table 5. Central Atomic^ and Molecular Hydrogen Densities 



Name 




n^2 from i^CO^ 


YiH^ from C^^C^ 


nffj/no 






(cm~^) 


(cm~^) 


(cm~^) 


from i^CO 


from C^^O 


L1544 


<3.0 


<3100 


6400 


4.8 X 10-^" 


2.3 X 10"^ 


B227 


2.0 


1000 


1200 


1.1 X lO--'^ 


8.3 X 10"^ 


L1574r 


3.9 


800 




2.4 X lO--'^ 




L 1574b 


5.9 


1100 


2400 


2.7 X 10-3 


1.2 X 10-3 


CB45 


6.0 


1700 


1900 


1.8 X 10-=^ 


1.6 X 10-3 



^Cold HI seen in absorption (HINSA) 

^Prom H2 column density derived from N(^3qo) given in Table 3, column 6 
^Prom H2 column density derived from N(C^^O) given in Table 3, column 7 
'^Total proton density in molecular and atomic hydrogen 
*^ Taking the source dimensions equal to upper limits measured 
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Table 6. Cloud Masses and Energies 



Name 


^vfwhm ^ 


^0^ 


M 


2T 


-U 




(km s~^) 


(10^ cm s-i) 


(1035 g) 


(10^5 erg) 


(1045 erg) 


LI 544 


0.64 


4.6 


4.8 


1.5 


5.7^ 


B227 


1.08 


7.9 


6.7 


4.8 


6.6 


L1574r 


1.42 


10.4 


1.9 


2.3 


0.76 


L 1574b 


1.43 


10.5 


7.6 


9.1 


8.5 


CB45 


1.50 


11.0 


13. 


17. 


25. 



^Full width to half maximum line width from ^3C0 

''3 dimensional velocity dispersion defined by (Tq = ['^{—^2^ — 7^)]°'^) 
where the line width and mass are those for ^3C0 

'^Based on Gaussian density distribution and central densities and dime- 
sions given in Table 4 and including a correction for helium 

'^Calculated using Gaussian density distribution; see Appendix B 

^'Assuming that wc have measured total mass, and that actual dimension 
is defined by our lower limit 



